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Abstract —This work explores the consequences of introducing 
a piezoelectric gate barrier in a normal field-effect transistor. 
Because of the positive feedback of strain and piezoelectric 
charge, internal charge amplification occurs in such an elec¬ 
tromechanical capacitor resulting in a negative capacitance. The 
first consequence of this amplification is a boost in the on- 
current of the transistor. As a second consequence, employing 
the Lagrangian method, we find that by using the negative 
capacitance of a highly compliant piezoelectric barrier, one can 
potentially reduce the subthreshold slope of a transistor below the 
room temperature Boltzmann limit of 60 mV/decade. However, 
this may come at the cost of hysteretic behavior in the transfer 
characteristics. 

Index Terms —Electrostriction, Electromechanical capacitor, 
Piezoelectric barrier, Negative capacitance, PiezoFET, Subthresh¬ 
old slope 

I. Introduction 

CALING of the size of field-effect transistors (FETs) has 
improved their performance and integration densities in 
integrated circuits for over two decades. Most conventional 
transistors make use of a passive insulating barrier layer 
between the gate metal and the semiconductor channel to 
modulate the density of conduction channel electrons or holes. 
Because the intrinsic properties of a passive gate barrier do 
not change with the applied voltage, they impose certain 
fundamental limitations on the resulting device performance. 

One such limitation is the subthreshold slope, i.e. the gate 
voltage required to change the drain current by an order of 
magnitude m, el given by 55 = m x 60 mV/decade at 
room temperature 0, (4|. Here, m = 1 + C sc /Ci ns is the 
‘body factor’, C sc is the semiconductor channel capacitance, 
and Ci n s is the gate insulator capacitance. In a traditional FET 
switch with a passive gate dielectric such as Si 02 , Ci ns > 0 
and thus m > 1, which leads to SS > 60 mV/decade g). This 
result, combined with circuit requirements for the on current 
I on and the on/off ratio I on /Ioff establish a minimum supply 
voltage Vdd, which does not scale in direct proportion with 
feature size ID, El, 0- Scaling of Vdd has hit a roadblock, 
giving rise to heat generation associated with the large power 
dissipation density in ICs 0, El, 0, 0, since the dissipated 
power is proportional to the square of the voltage, Pdiss oc V‘ r jd 
0, 0- Many ideas based on alternate transport mechanisms 
in the semiconductor channel, such as interband tunneling, or 
impact ionization are being explored to lower Vdd- 

An interesting alternative is to replace the passive gate 
barrier with an active one. A first proposal of an active 
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ferroelectric insulator 0 predicts internal voltage gain: the 
voltage across the gate insulator layer is larger than the applied 
external gate voltage. The origin of internal voltage gain is 
the collective alignment of the microscopic electric dipoles 
in the ferroelectric layer in response to the external electric 
field produced by the gate voltage. The alignment of dipoles 
generates a voltage of its own, thus amplifying the voltage 
that makes it to the semiconductor channel. Under appropriate 
bias conditions 0 , the insulator capacitance provided by the 
ferroelectric is mathematically negative ( Ci ns < 0), causing 
to = 1 + C ac /Ci ns < 1 and SS < 60 mV/decade. Such 
an active-gate FET then will require a lower gate voltage to 
create the same charge as a conventional FET with passive 
gate dielectrics 0 , thereby facilitating device scaling. 

In this paper, we explore the device consequences of us¬ 
ing a piezoelectric insulator as the active gate barrier in a 
transistor instead of the ferroelectric barrier. Piezoelectric gate 
barriers are at the heart of commercially available Ill-nitride 
heterostructure transistors 0, 0- We first consider an active 
compliant piezoelectric layer as the insulator in a parallel plate 
capacitor. We find that this simple electromechanical capacitor 
system exhibits a remarkably rich range of behavior. We show 
that negative capacitance emerges as a natural response to 
applied voltage. In this regime of negative capacitance, we 
show that we obtain a higher charge than in a corresponding 
capacitor with a passive dielectric. Non-trivial capacitance- 
voltage behavior in such capacitors have also been reported 
experimentally 01, ED- Next, we port the parallel-plate 
electromechanical capacitor to the gate capacitor of a FET. 
We show how this piezoelectric gate stack enables a higher 
on-current than in a transistor with a passive dielectric due 
to internal charge amplification. Finally, building upon our 
earlier proposals ED, ED, ED, we discuss the possibility of 
using the negative capacitance regime of a highly compliant 
piezoelectric barrier to obtain sub-60 mV/decade switching in 
a transistor. 

II. Electromechanical capacitor 

We begin by discussing the piezoelectric parallel-plate ca¬ 
pacitor. Consider the parallel plate capacitor of area A shown 
in Fig. |T](a). The equilibrium thickness to of the piezoelectric 
insulator layer sandwiched between the metal plates changes to 
t = to — 5 when a voltage V is applied on the plates, as shown 
in Fig.[H(b). The strain is defined as s = 6/to. The equal and 
opposite sheet charges o m that develop on the metal plates 
set up an attractive force between them, which strains the 
insulator. This effect, called electrostriction, is the electric-field 
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to the metal charge lfl5l . lH 6 l via P = F/A = cr^/ed- To 
linear order, the pressure depends on the strain via the stiffness 
coefficient P = C 33 S, where C 33 is in units of N/m 2 , or 
Pascals. Thus, we obtain the strain as a function of the metal 
charge: s = e ( iC:>, 3 . Substituting in Eq. |T] and rearranging, 

we have the desired relation between the metal charge in 
response to an applied voltage: 


Fig. 1. (a) Schematic cross section of a parallel-plate electromechanical 

capacitor with piezoelectric barrier layer of thickness to at V = 0 Y, b) The 
layer thickness shrinks to to — 5 when voltage V is applied. Sheet charge 
distribution p(z ) with ±<7 m on the metal plates and surface charges =b<r s on 
the piezoelectric. 


induced reduction of the thickness of a material; it occurs in all 
insulators, whether or not the layer is piezoelectric. However 
if the insulator is piezoelectric, the strain amplifies the surface 
charge of the insulator. This mechanism sets up a positive 
feedback between the thickness and the electric field, and 
is responsible for the appearance of negative capacitance. To 
find the capacitance in the presence of such electromechanical 
coupling, one must first find the net metal charge <r m as a 
function of the external (battery) voltage, and then take its 
derivative. This requires us to identify the surface charges 
a s that develop at the surface of the insulator. The resulting 
electric field profile is constant, equal to E = V/t, and the 
voltage drops linearly across the insulator. 

Maxwell’s boundary conditions across the metal-insulator 
interface requires the normal components of the displacement 
vector to obey Dd — D m = a m . Dd is the displacement 
field in the dielectric related to the surface charges a s by 
D d = e 0 E + a s , where a s = ( e d - e 0 )E + e 33 s + a sp or 
D d = e d E + e 33 s + a sp . Here e d = e 0 (l + Xd) is the net 
dielectric constant of the piezoelectric layer, and Xd is its 
electric susceptibility. The electric field is E = V/t, where 
t = fo(l ~ s) is the thickness of the strained insulator layer. 
We explicitly allow for both piezoelectric and spontaneous po¬ 
larization for an active dielectric material. The strain-induced 
piezoelectric contribution to the charge (to linear order) is 
e 3 3 S, where e 3 3 is the piezoelectric coefficient in units of C/m 2 
and s = 6 /to is the strain along the field. The charge due to 
spontaneous polarization is a sp , also in units of C/m 2 . Inside 
the metal, D m = 0. Therefore, we obtain the relation 
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This relation illustrates how the strain s explicitly enters 
the electrostatic relation between the metal charge and the 
voltage across the capacitor. If one neglects the spontaneous 
polarization ( <r sp —> 0 ), piezoelectric effect (e 33 — > 0 ) and 
strain (s — > 0), we get a m = CqV, (with Co = e d /to), 
the standard textbook formula of a parallel plate capacitor. 
However, we note that one can turn off the spontaneous and 
piezoelectric polarization by choice of material, and yet the 
factor (1 — s) in the denominator will persist: this is the 
electrostriction term. 

The mechanical pressure P experienced by the insulator 
is the electrical force F per unit area A. It is thus related 
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The right hand side is a fourth order polynomial in a m , 
and captures the electromechanical coupling physics. Let us 
explore its consequences. The sheet charge on the metal 
n m = a m /q from Eq. [2]is plotted as a function of the applied 
voltage V in Fig. [2] for different sets of material parameters. 
For example, e 33 = 3.1 C/m 2 , e d = 15eo correspond to the 
piezoelectric material Sc x Ali_ x N fTTI . JT 8 ). The value of C 33 
is allowed to vary arbitrarily in order to investigate the range 
of behavior of the piezoelectric capacitor. We also assume 
( 7 sp = 0. A non-zero value of cr sp merely causes a horizontal 
shift of the a m — V curve (see Supporting document). The 
physics of the piezoelectric capacitor with cr sp = 0 becomes 
apparent by factoring Eq. [2] into 
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where qn m = a m , qn v = e d C3 3 /e33, and qn v = \Ze d C 3 3. 
Setting V = 0 in Eq. [3] we obtain four real roots n m 0 = 0, 
+ 1 irj, —n v and Vn^. For a rigid (C 33 —y 00 ), non-piezoelectric 
(e 33 = 0 ) insulator, n^^n-n 00 , whereupon we recover 
<r m = qn m = CqV, and the metal charge is a linear function 
of voltage as shown in the green line in Fig. [3 a). 

On the other hand, for a compliant non-piezoelectric in¬ 
sulator, C 33 > 0, and Eq. [3] reduces to a cubic equation 
with roots 0, ±n v at V = 0. This is in fact a prototypical 
description of a nano-electromechanical switch lfl9l . The two 
additional roots ±n v make the dependence of a m on V 
nonlinear with two additional zero crossings. Multiple zero 
crossings of the qn m — V curve mathematically guarantees 
that there must be regions of negative slope d(qn m )/dV < 0. 
This is shown in red in the flipped S-shaped curve of Fig. 
|2|a), where C 33 = 1 GPa is assumed. In these regions, the 
electromechanical capacitor has a negative capacitance. 

However, the negative capacitance corresponds to very high 
values of charge density(> 10 14 cm -2 ) and strain s (> 0.3), as 
shown in Fig. |2a). Though rapid progress is being made in the 
solid-state electrostatic gating of ever increasing carrier densi¬ 
ties in semiconductors l 20 l . methods to reduce the charge and 
strain are desirable. Now consider the piezoelectric insulator 
where C33 > 0 and e 3 3 > 0. Notice that n v is independent of 
e 3 3 . The root n n depends on 633 and its location determines 
the shape of the qn m — V curve. If e 3 3 > \/c d C 33, then 
0 < n-n < n v and negative capacitance appears in the two 
charge segments [n\ = —n v ,n,2] and [ 713 , 714 ] shown in red 
in Fig. [2] It is important to realize that piezoelectricity lowers 
both the charge density (~ 10 14 cm -2 ) [inset of Fig. [3a)] 
and strain < 0.01 [Fig. [3jb)] at which negative capacitance 
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Fig. 2. a) Charge-voltage (qn m —V) characteristic of the electromechanical capacitor. Various charge states such as the positive capacitance segments [712, ?':+ 
+4, 77.5] where the slopes Cpe = d(qnm)/dV > 0 are positive, and negative capacitance segments \ri \, rr2], [713,714] where Cpe = d(qn rn )/dV < 0 are 
shown, b) The characteristics of a piezoelectric capacitor with a lower stiffness and more compliant barrier with C33 = 0.01 GPa makes negative capacitance 
accessible at a lower charge ~ 10 11 cm -2 , as shown in the inset. 



Fig. 3. Strain as a function of voltage in a) a non-piezoelectric insulator 
layer with C33 = 1 GPa, 633 = 0 C/m 2 , to = 0.5 nm, and b) in a 
piezoelectric insulator layer with C33 = 1 GPa, 633 =3.1 C/m 2 . Strain s < 1 
is physically accessible in solid state, where the remaining layer thickness, 
/,Q (I — s) > 0. Red (black) represents strain corresponding to the negative 
(positive) capacitance charge states. 


appears, compared to electrostriction alone. For vanishingly 
small voltages around zero, the piezoelectric capacitor behaves 
exactly like a parallel plate capacitor - a straight line. But 
the additional benefit of the above coupling is the increased 
charge density compared to a passive dielectric due to the 
piezoelectric amplification - this effect will boost the on- 
state current in a transistor. From Eq. Q] the piezoelectric 
amplification is n m — ^- « + to leading order. Finally, 

if we use a highly compliant piezoelectric, for example with 
C 33 = 0.01 GPa, negative capacitance can be accessed at very 
low charge density ~ 10 11 cm -2 , as indicated in Fig. |2jb). 
These highly compliant piezoelectrics can potentially enable 
the design of transistors with steep sub-threshold behavior, 
but require new materials as will be described later. We also 
remark here that Pauli’s exclusion principle of solid matter 
and quantum compressibility restricts s < 1. Therefore for 
piezoelectric insulators, the metal charge will be restricted to 


—< n m < + 71 ,,. It may be possible to go beyond these 
restrictions (s > 1, shown as dashed lines in Fig.|2j in gaseous 
plasmas where charged ion plate ‘electrodes’ can pass through 
each other. But we do not pursue that line of analysis here, 
by restricting the discussion to solid metals and dielectrics. 

III. Transistor with a Piezoelectric Barrier 




Fig. 4. Schematic cross section of a transistor (“piezoFET”) with a 
piezoelectric gate barrier, semiconductor channel such as Si, GaN or 2D 
material, M0S2, and source and drain ohmic contacts. The gate capacitance 
circuit is a series combination of the piezoelectric capacitance Cpe and 
the semiconductor capacitance Csc- Here intrinsic gate voltage V' ga = 
Vgs — !,[ Rs, and intrinsic drain voltage Vg = V+ — /)/ ( R s + Rg) . where R a 
and //,, are the source and drain contact resistances. The energy band diagram 
is shown for the metal-piezoelectric-semiconductor stack of the transistor. 
i/j s = Vg S — V = ( E — /+;) l([ is the surface potential. 

We now explore how the presence of the piezoelectric 
capacitor in the gate of a transistor with a semiconductor chan¬ 
nel changes the traditional characteristics. The semiconductor 
channel could be formed of a gapped 2-dimensional crystal 
such as M 0 S 2 , or a 3-D crystal semiconductor such as Si or 
GaN. The semiconductor is characterized by the valley degen¬ 
eracy g v of the conduction (or valence) band. We assume the 
energy dispersion of each valley to be the same, characterized 
by an effective mass m* and spin degeneracy g s = 2. Carrier 
transport in the semiconductor channel is assumed to be 2- 
dimensional - which holds both for monolayer 2D crystals 
and in field-effect transistors made of 3D semiconductors. 
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Fig. 5. a) Gate capacitance C g versus V ga for transistors with piezoelectric (solid line) and dielectric (dashed line) insulators, b) Graphical load line analysis 
to obtain sheet carrier density n s for different V gs . Blue curves show the semiconductor charge for different V gs , green shows the metal charge in the case 
of a passive dielectric, whereas red(black) shows the metal charge in the negative(positive) capacitance regimes of the piezoelectric capacitor. Intersections 
ai and <12 of the above characteristics define the operating points of the system, c) Transfer curve depicts the drain current I d versus V gs at drain voltage 
V da = 0.5 V for GaN transistors with piezoelectric (solid line) and dielectric barriers (dashed line). 


where transport occurs in a quasi-2D electron/holes gas. The 
occupation of multiple 2D subbands can then be treated as 
individual 2D channels - we consider a single subband model. 

The semiconductor channel of length L and width W 
is assumed to be connected to very low-resistance ohmic 
contacts at the source and drain, as shown in Fig. 0] The 
energy-band diagram in Fig. 0] shows the potential barrier 
controlled by the voltage on the gate metal. Electrical charge 
neutrality requires <j m = qn s , where n s is the mobile 
carrier sheet density at the ‘top-of-the-barrier’ in the energy 
band diagram along the length of the channel. The energy 
band diagram from the metal to the semiconductor requires 
qcpB + qV — A E c + (E Fs ~ E c ) = qV gs . By suitable choice 
of materials, we assume that qifip = A E c \ if this is not the 
case, the difference can be absorbed in a shift of threshold 
voltage. When no drain voltage is applied, carriers in the 
semiconductor are in thermal equilibrium with the source and 
drain reservoirs, which for a parabolic 2D bandstructure means 
qn s = C sc V th In (1 + exy>[(E Fs - E c )/kT\), or E Fs - E c = 
kT In (exp[qn s /C sc V t h] - 1), where C sc = q 2 g s g v m*/2 ttH 2 
is the density of states semiconductor capacitance, and the 
thermal voltage Vth = kT / q. From the energy band diagram 
in Fig. 0] the relation between the applied gate voltage V gs 
and the voltage drop V across the piezoelectric insulator is 
qV' s =qV+ ( E Fs - E C ). Here, (E Fa - E c )/q = (V' s - 
V ) = ip s is the surface potential. Using the carrier density 
expression and Eq. 0] the gate-induced charge qn s in the 
semiconductor channel is self-consistently calculated. Finally, 
using this new dependence of charge on the voltages and the 
piezoelectric coefficients, the current-voltage characteristics of 
the piezoFET are obtained from the ballistic transport model 
ED incorporating the quantum contact resistances of 0.026 
kO.//m ll22l at the source and drain ends. 

Fig- 0 shows the gate capacitance C g = d(qn s )/dV gs , and 
device characteristics (I — V) of a ballistic piezoFET with a 


GaN channel (to* = 0.2 mo, g v = 1 E3). Fig. 01a) shows 
that a higher gate C g is obtained in the piezoFET (solid line), 
as compared to a FET with a passive gate (dashed line). The 
higher C g is due to the negative capacitance resulting from 
piezoelectric charge amplification: CpeC sc /(Cpe + C sc ) > 
C 0 C sc /(C 0 + C sc ) when Cpe < 0. Fig.01b) depicts the solu¬ 
tion of the piezoelectric and semiconductor charge equations 
graphically, following the load-line approach (see 1241). The 
blue lines depict charge in the semiconductor channel, and the 
green, black, and red lines depict the charge drawn into the 
metal from the battery. They must be equal to maintain global 
charge neutrality, meaning the locus of intersections are the 
operating points of the device. The green line is the charge 
on the metal for a traditional passive gate dielectric, and the 
red/black lines for a piezoelectric gate. When the transistor is 
on (V' ~ + 0.3 V), an increase in the charge at point <22 in 
Fig. 01 b) is seen for the piezoelectric compared to point aq 
for a passive dielectric. This increased charge boosts the on- 
current as depicted in Fig. 01c), consequently improving the 
Ion/Ioff ratio. This sort of piezoelectric amplification is an 
interesting method to boost the on-current in any transistor. 
Since much of the high-performance characteristics such as 
gain and cutoff frequencies depend on I on , corresponding 
boosts can be expected in these parameters. This may be 
specially useful for boosting the current in FETs made of 
relatively low mobility channel materials. Note however in 
Fig. 01c) that this device still has a SS of 60 mV/decade. This 
is because the negative capacitance regime is only accessible 
for charge densities > 1.5 x 10 13 cm -2 : at this high level 
of charge, the transistor is in its on-state, rather than in the 
sub-threshold regime. 

Because the charge-voltage characteristic of the piezoelec¬ 
tric capacitor is highly non-linear, it can have multiple inter¬ 
sections with the semiconductor load line. Ref. l l25l develops 
a systematic procedure to understand such non-linear systems 
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Fig. 6 . a) Load line analysis showing multiple intersections of the piezoelectric and semiconductor characteristics for different V ' s , b) Free-energy landscape 
of the piezoelectric-semiconductor stack at various V' s . Blue and black dots show stable operating points. 



Fig. 7. a) Id versus V gs curves at Vd s = 0.1 V for a GaN channel piezoFET with a compliant piezoelectric with C 33 = 0.01 GPa. A boost in the on-current 
and sub-60 mV/dec SS (inset) are obtained as compared to a passive gate dielectric, b) Load line characteristics to explain the hysteresis with gate bias 
voltages Vg S , c) The calculated hysteresis in the transfer curve Id versus V gs for forward and reverse sweeps is shown for a GaN channel piezoFET. 


based on the Euler-Lagrange equations of motion (see support¬ 
ing document for details). For this analysis, we define a free- 
energy G in units of J/m 2 for the piezoelectric-semiconductor 
stack: G(a m , V' s ) = f Vda m +f ip s da m -^mV' s , where V is 
the voltage drop across the gate insulator and i/) s is the surface 
potential. Minima in this free-energy landscape correspond to 
stable charge solutions of the non-linear system. If there are 
multiple minima, the actual solution cr m = qn s depends on 
the previous state, or the history of the system. 

For example. Fig. [6|a) shows the load lines and the corre¬ 
sponding evolution of the free-energy landscape for different 
Vg S are shown in Fig (bj'b). The shape of the energy landscape 
changes with the applied voltage. There are two energy 
minima in the range —0.35 < V' js < 0.35 V, and a single 
minimum otherwise. Let us assume that there is no charge to 
begin with on the capacitor, and ramp the gate from a negative 


to a positive voltage. Until around V' s = + 0.35 V, the system 
remains in the minimum corresponding to the lower charge 
state (~ 2.3 x 10 13 /cm 2 ) shown as a blue dot in Fig. [6ja) and 
the inset of Fig. [6jb). But when Vg S > 0.35 V it is driven into 
the higher charge state shown as a black dot. Thus, provided 
Vg S < 0.35 V, the transistor displays no hysteresis in its I — V 
characteristics. 

It is pertinent here to note an important difference in 
the nature of negative capacitance of the piezoelectric and 
ferroelectric insulators. The ferroelectric capacitor possesses 
negative capacitance at zero charge, whereas the capacitance 
of the piezoelectric capacitor is positive at zero charge. This 
property of the ferroelectric capacitor is exploited in achieving 
SS < 60 mV/decade, since the semiconductor load line can 
intersect the negative capacitance regime of the ferroelectric 
characteristic at the very low charge densities corresponding to 
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subthreshold operation of the transistor. Can a similar negative 
capacitance be obtained in the SS regime (Vg S < 0 V) using 
piezoelectric gates? We explore this by tuning the piezoelectric 
material properties. 

We find that if a lower stiffness, highly compliant piezo¬ 
electric barrier with C 33 ~ 0.01 GPa is used, it can enable 
the reduction of the subthreshold slope below 60 mV/decade 
and also boost the on-current. This is shown in Fig. [7] a). 
Here negative capacitance is accessed in the subthreshold 
region, shown by the operating point ai in the load line 
characteristics at V' ga = -0.05 V shown in Fig. |7Jb). The 
on-state operation of this transistor corresponds to the higher 
charge state determined by the operating point a 2 in the 
load line characteristics at V' =0.1 V. However, this also 
results in hysteresis in the transistor characteristics with V gs 
sweep, as shown in the Id — V gs characteristics in Fig. |7Jc) 
which is calculated using the Lagrangian method. Hysteresis is 
undesirable in purely switching applications, but desirable for 
memory. Further, the strain in the higher charge states < 12 , «3 is 
very close to 100 %, which is not feasible in realistic materials. 
If suitable new piezoelectric materials with ultra-low C 33 and 
high 633 could be developed (see Supporting Information 
for various piezoelectrics with different C 33 and 633 ), sub- 
60 mV/decade switching can be achieved with hysteresis 
with suitable choice of semiconductors. Investigation of other 
transistor designs incorporating the piezoelectric barrier, such 
as the quantum metal transistor ||26l to eliminate the hysteresis 
and reduce strain could be the focus of future work. 

IV. Conclusion 

We also emphasize that we have assumed linear piezoelec¬ 
tric parameters in this work to keep the model simple and 
yet capture the new physics. The non-linear material response 
needs to be explored in future. To conclude, the behavior 
of transistor switches using active piezoelectric gate barriers 
was explored. Because of electrostriction and piezoelectricity, 
negative capacitance is predicted to appear in a piezoelectric 
capacitor. Using this negative capacitance and a ballistic trans¬ 
port model, we predict that compliant piezoelectric barriers 
can boost the gate capacitance and increase the on-currents 
of transistors. Also, steep switching with sub-60 mV/decade 
subthreshold slope is predicted when the negative capacitance 
of the piezoelectric barrier is accessed in the off-state operation 
of the transistor, and this steep behavior is predicted to be 
assisted by hysteresis based on the Lagrangian method of 
stability of the transistor system. 
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SI. (Jm-V RELATION OF PIEZOELECTRIC CAPACITOR WITH <7. 


sp 


Following the same notation as the main text, the 4 th order charge versus voltage a m — V relation of an electromechanical 
capacitor is 

(a 3 p — 633 ) 2 1 3 , e 33 4 


CqV — Um Us 


£d.C 33 


-OVn — 


tdC 33 


{^dC 33 ) 


2 m- 


(Si) 


The right side is a fourth order polynomial in u m , and captures the electromechanical coupling physics. Writing u m = qrir, 
where n m is the sheet charge density on the metal and a sp = qn sp ^ 0, the polynomial factorizes to 


C 0 V = qn sp ( 1 


- 1)(1 - —)(1 + —), 

Tlct Tla ^77 ^77 


(S2) 


where n+ = (n , r + \/n\- 4n sp n n )/2, n a = (n n - s/n*- 4n sp riT r )/2, qn n = e d C 33 /e 33 , and qn v = \Ze d C 33 . The 
four roots are characteristic sheet densities determined uniquely by the electromechanical coefficients and the spontaneous 
polarization of the dielectric material. If > 4 n sp which is met if C 33 > 4n sp e 33 /e d , then all four roots are real. The effect 
of spontaneous polarization a sp ^ 0 C/m 2 is a voltage offset, which leads to left and right shifts of the qn m — V characteristics 
(black and red curves) with respect to the characteristic (blue curve) with u sp = 0 C/m 2 , as shown in Fig. ED These shifts 
will move the threshold voltages of a corresponding transistor, and will also locally change the slopes of the charge-voltage 
characteristics. 



Fig. SI. Charge-voltage (qn m — V ) characteristics of an electromechanical capacitor with piezoelectric barrier. Color lines (black, blue and red) show the 

characteristics of capacitors for a 3p = 0.2 C/m 2 , 0 C/m 2 , and -0.2 C/m 2 . a ap 7 ^ 0 C/m 2 leads to horizontal left and right shifts of qn m — V curve. 


S2. Ballistic FET I-V model 

When no drain voltage is applied, carriers in the semiconductor are in thermal equilibrium with the source and drain 
reservoirs, which for a parabolic 2D bandstructure means qn s = C sc Vth In (1 + exp[(E/F S — E c )/kT)). From the energy band 
diagram of metal-piezoelectric-semiconductor stack of a transistor (Fig. IS2I) . the voltage division qV gs = qV + (Ep a — E c ) 
translates to the dimensionless equation, 

V gs qF(n s ) qn s 

e^th = e c o v th (e.v q v t h — l) ; (S3) 

where we have defined F(n s ) = n sp {l — -^)(- 2 =r —1)(1— ^■)(1+^ 1 ). This equation must be solved to find the semiconductor 
charge n s at a gate voltage V gs . If the dielectric is piezoelectric but does not have spontaneous polarization, one must replace 
F[n s ) by F pz (n s ) = n s (l — —)(1 + is -)(l — ■ 2t ) to find the semiconductor charge in response to the gate voltage. This is 
the major change to a standard Natori-type ballistic FET model ID brought about by the piezoelectric gate barrier. 

When a drain voltage V ds is applied, the carrier distribution in the ‘top-of-the-barrier’ point x rna:r in the energy band diagram 
is split in two. In the ballistic limit of transport, the right-going carriers are in equilibrium with the source reservoir of Fermi 
energy Ep s , whereas the left-going carriers are in equilibrium with the drain reservoir Ep d , and they are out of equilibrium 
by Ep s — Ep d = qV ds . Note that Ep s ^ Ep s \ the application of a drain bias causes a rearrangement of the carrier distribution 
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Fig. S2. Schematic cross section of a transistor (“piezoFET”) with a piezoelectric gate barrier, semiconductor channel such as Si, GaN or 2D material, 
M 0 S 2 and source and drain ohmic contacts. Energy band diagram for the metal-piezoelectric-semiconductor stack of the transistor. For transistor operation, 
we use ballistic transport model to calculate the transistor characteristics. The mobile sheet carrier density at the 'top-of-the-barrier' in energy band diagram 
is controlled by V gs through the piezoelectric gate barrier, and source and drain Fermi levels Ep s & Ep d are separated by . Dashed circle shows the 
/r- space for carrier distributions at Va s = 0 and grey circle shows the fc-space for carrier distributions at applied V d i.e. Ep a — Ep d = qV d , v g (k) is the 
group velocity of carriers. Here, V dg = Vds ~ Id(Rs + Rd)> where R s and R d are the source and drain contact resistances. 


in the semiconductor channel. However, with good electrostatic design, one can ensure that the net carrier density at x max is 
the same as when V ds = 0. The carrier distribution in the k —space is depicted in Fig. IS21 the dashed circle is the distribution 
for V ds = 0, and the gray half-circles are the result of application of a drain bias, both in the T —> 0 K limit. Defining 
rj s = ( Ep s — E c )/kT, Vd = V ds /kT, we find that for maintaining the same carrier density, one must meet the condition 

qn s = \C q Vth ln(l + e ris )(l + e r,B ~ Vd ), which yields 


Vs 



2gn s 

_|_ 5^)2 _|_ 4 e «d( e c,v ik 


!) — (! + e Vd )\ — In [2] 


(S4) 



Fig. S3. Output characteristics depicts the drain current Id versus Vd s at gate bias voltage V gs = 0.2 V with a step of -0.05 V for GaN piezoFETs with 
piezoelectric (solid line) and dielectric barriers (dashed line). A boost in the on-current is obtained for transistor with a compliant piezoelectric barrier than a 
dielectric one. 


Note that this is of a similar form as the seminal result on ballistic transistors by Natori et al (TJ, but the factor n s inside 
the square root must be obtained from Eq. [S3] to account for the electromechanical coupling self-consistently. By summing 
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Fig. S 4 . Piezoelectric Materials with C'3,3 and 633 coefficients. Table for different C33 and 633 coefficients of various piezoelectrics with corresponding 
references. C33 and 633 values are obtained by literature review from listed references. 633 is extracted from a matrix relation of [e] = [C] * [d] where [e] 
(unit: C/m 2 ) is the piezoelectric coefficient matrix, [C] (unit: N/m 2 ) is the elastic stiffness matrix, [d] (unit: C/N) is the piezoelectric coupling matrix. These 
matrices have sparsity pattern depends on the symmetry of the crystal. However, here in order to get an estimation of the magnitude of e33, we use a simpler 
relation 633 (C/m 2 ) = C33 (N/m 2 ) * d 33 (C/N), in cases where 633 is not quoted in references. The arrow shows the desired piezoelectric materials to be 
used for gate barriers of steep transistors. 


over the group velocities of the k— states, the net current per unit width of the ballistic piezoFET is then given by the same 
expression as in Natori (TJ: 

J= Jo (>1 fo) - Fi ( 77 , - %)) (S5) 

_ 3 

where Jo = ggsg " ’ anc * * s t ' le Fermi-Dirac integral of order 1/2. 

The above expressions provide the complete electrical characteristics of a ballistic piezoFET at any temperature in a compact 
model. We can obtain n s at any gate voltage V gs fully accounting for the electromechanical coupling by solving Eq. [S3] We 
then find r/ s at any given V gs and V,i s using Eq. [S4] and finally find the drain current per unit width using Eq. [S5] Fig. [S3] 
shows the output characteristics depicting /,/ — V r / S at different V gs for GaN channel transistors with a compliant piezoelectric 
(solid lines) and passive dielectric (dashed lines) barriers. 

Fig. [S4] describes different piezoelectric materials with C 33 and 633 coefficients obtained by literature review from Refs. 0 , 
a, a, a, a, a, 0, 0. Compliant piezoelectric materials with parameter space such as lower C 33 (higher compliance), 
and higher e 33 (higher piezoelectricity), shown by the arrow in Fig. [S4] are useful for gate barriers of proposed steep transistors. 


S3. Energy landscape for piezoelectric-semiconductor stack: Computational Details 


A. Definition of Free-Energy 

Consider the circuit in Fig. [S5] showing a piezoelectric-semiconductor stack connected to a voltage source via a resistor R 
(units: flcm 2 ). Both the piezoelectric insulator and semiconductor have non-linear charge ( a m ) - voltage ( V) characteristics, 
denoted say by V = V/ rls = fi(a m ) and s = / 2 (cr m ) respectively. A systematic method of studying the behavior of circuits 
with non-linear elements uses the Euler-Lagrange equations GO): 


■£(^ra 1 Gm) 
U{(Jm) 

r<7m 

+ / f2((Tm)d<Tm 
JO 

dam ( dC \ dC 


/ OL, \ 

V d<7 rn J 




T' -U 

/ fi(vm)da n 


= V n , 


(S 6 ) 


(S7) 


dt 


— f?cr, 
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System 



Fig. S5. Schematic of a system shows a piezoelectric-semiconductor stack, represented by a series combination of piezoelectric capacitance Cpe and 
semiconductor capacitance C sc connected to V gs via a resistor R. ip B is the channel surface potential at the insulator-semiconductor interface. 


where £(cr m ,cr m ) is the Lagrangian, T'(& m ) is the magnetic co-energy (in any inductors) and U(a m ) is the potential energy 
in the capacitors in the circuit. For the circuit in Fig. IS5I T' = 0, and we have 

h{(Jm) + f2{cJ m ) = V' gs -Rdm, (S9) 

which under equilibrium (& m = 0) gives V' gs = /i(cr mo ) + f2(o’mo). To determine whether a mo is a point of stable/unstable 
equilibrium, we write <j m = &m 0 +<5cr mo for a small perturbation 6 a mo , expand the Lagrangian as a Taylor series (upto second 
order) about a mo , and analyze whether the perturbation grows or decays with time: 

S am = _ 2 /i(^o)+/^(^o ) S(Tm = Sa m /r (S10) 

Scr m {t) = 6a mo exp(-t/T). (SI 1) 

where f'(cr m ) is the derivative of f(a m ) w.r.t a m . The system is stable to perturbations if 6 <j m (t ) —► 0 as t —> 0, i.e. if 
fi( a m 0 ) + /^mj > 0. The above discussion motivates the definition of the free-energy G(o m ): 

f‘&m 

G(cr m , Vg S ) = / [fl(&m) + /2(^m)] dam ~ °mVg S - (S12) 

Jo 

From the carrier density expression in a semiconductor channel, ip s can be expressed as ip s = f'li&m) = 
Vth In (e c * cV th — 1^. Using the expressions of ijj s and Vi ns = in Eq. IS 1 21 the free energy G (unit: J/m 2 ) is found to 

be 


G — 0 : 5 ( 7 ^ + 0 : 4 ( 7 ^ + 

+V th j In ^e c ^ v t-h - da rn - (J rn V gs , (S13) 


where «5 = 5 c 0 (^p 33 )a , 0:4 = 4 Coe 1 d ig 33 , «3 = 3 Coejc 33 ’ 012 = 2 are mater ial and geometric constants of the problem. 
Note the internal energy is a non-linear function of the sheet charge a m , and is a linear function of applied gate bias voltage 
V gs . Then, all points of equilibrium satisfy dG/da m = 0, and further, points of stable equilibrium satisfy d 2 G/daf rl > 0. 


B. Self-consistent solution 

We make use of the Euler-Lagrange equation to calculate the I —V characteristics of the PiezoFET. The basic structure of 
the self-consistent algorithm involves an outer loop that determines the voltage drop across the intrinsic FET (i.e eliminating 
the voltage drops across the contacts). Inside this loop, we require the charge n s for the intrinsic gate voltage. To do this, 
we solve the Euler-Lagrange equation (using the implicit Euler method for about 5r) to determine a good guess of charge 
n s guess at the present voltage, starting from the value of n s 0 obtained at a previous voltage. n Stguess is then used in a fixed 
point iteration scheme to determine n s , from which the current is finally calculated. We have not encountered prior use of this 
Euler-Lagrange method for the calculation of hysteretic characteristics of electron device systems in an extensive literature 
search, and intend to publish the detailed procedure in a follow-up report. 
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